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ABSTRACT 

The correlation between the peak photon energy of the internal spectrum E Pi i and 
isotropic equivalent radiated energy Ei so (the Amati relation) is explored in a scalar 
field model of dark energy. Using an updated data set of 109 high rcdshift GRBs, we 
show that the correlation parameters only weakly depend on the cosmological model. 
Once the parameters of Amati relation have been determined we use this relation to 
construct a fiducial GRBs Hubble diagram that extends up to redshifts ~ 8. Moreover 
we apply a local regression technique to estimate, in a model independent way, the 
distance modulus from the recently updated Union SNIa sample, containing 557 SNIa 
spanning the redshift range of 0.015 < z < 1.55. The derived calibration parameters 
are used to construct an updated GRBs Hubble diagram, which we call the calibrated 
GRBs HD. We also compare the fiducial and calibrated GRBs HDs, which turned 
out to be fully statistically consistent, thus indicating that they are not affected by 
any systematic bias induced by the different calibration procedures. This means that 
the high redshift GRBs can be used to test different models of dark energy settling 
the circularity problem. Furthermore, we investigate possible evolutionary effects that 
might have important influence on our results. Our analysis indicates that the presently 
available GRBs datasets do not show statistically unambiguous evolutionary effect 
with the cosmological redshift. Finally we propose another approach to calibrate the 
GRB relations, by using an approximate luminosity distance relation, which holds in 
any cosmological model. We use this calibration of the Amati relation to construct 
an empirical approximate HD, which we compare with the calibrated GRBs HD. We 
finally investigate the implications of this approach for the high redshift cosmography. 

Key words: Gamma Rays : bursts - Cosmology : distance scale - Cosmology : cos- 
mological parameters 



1 INTRODUCTION 

At the end of the '90s observations of high redshift supernovae of type la (SNIa) revealed that the universe is now expanding 
at an accelerated rate. This surprising result has been independently confirmed by observations of small scale temperature 
anisotropies of the cosmic microwave background radiation (CMB) (Riess et al. 2007; Asticr ct al. 2006; Kowalski ct al. 2008 
|Spergel et al. 20 07). It is usually assumed that the observed accelerated expansion is caused by a so called dark energy, with 
unusual properties. The pressure of dark energy pde is negative and it is related to the positive energy density of dark energy 
tde by pde = w ede where the proportionality coefficient w < 0. According to the present day estimates, about 75% of matter- 
energy in the universe is in the form of dark energy, so that now the dark energy is the dominating component in the universe. 
The nature of dark energy is not known. Proposed so far models of dark energy can be divided, at least, into three groups: a) a 
non zero cosmological constant, in this case w — — 1, b) a potential energy of some not yet discovered scalar field, or c) effects 
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connected with non homogeneous distribution of matter and averaging procedures. In the last two possibilities, in general, w 
is not constant and it depends on the redshift z. Observations of type la supernovae and small scale anisotropies of the cosmic 
microwave background radiation are consistent with the assumption that the observed accelerated expansion is due to the 
non zero cosmological constant. However, so far the type la supernovae have been observed only at redshifts z < 2, while in 
order to test if w is changing with redshift it is necessary to use more distant objects. New possibilities opened up when the 
Gamma Ray Bursts have been discovered at higher redshifts, the present record is at z = 8.26 ijGreiner et al. 2009ft . GRBs arc 
however enigmatic objects. First of all the mechanism that is responsible for releasing the incredible amounts of energy that 
a typical GRB emits is not yet known (see for instance Meszaros 2006 for a recent review). It is also not yet definitely known 
if the energy is emitted isotropically or is beamed. Despite of these difficulties GRBs are promising objects that can be used 
to study the expansion rate of the universe at high redshifts ( |Bradley 2003] ISchaefer 20031 IDai et al. 20041 IBloom et a~ 2003 
IFirmani at al. 20051 ISchaefer 20071 ILi et al. 20081 lAmati et al. 20081 ITsutsui et al. 2009f> . Using the observed spectrum and 
light curves it is possible to derive additional parameters, for example the peak photon energy E Pt i, at which the burst is 
the brightest, and the variability parameter V which measures the smoothness of the light curve (for definitions of these and 
other parameters mentioned below, see Schaefer (2007)). From observations of the afterglow it is possible to derive another 
set of parameters, redshift is the most important, and also the jet opening angle Qjet, inferred from the achromatic break in 
the light curve of the afterglow, the time lag r; a9 , which measures the time offset between high and low energy GRB photons 
arriving at the detector, and trt - the shortest time over which the light curve increases by half of the peak flux of the burst. 
Moreover, even though the most important parameters - the intrinsic luminosity L and the total (isotropic) radiated energy 
-Eiso - are not directly observable, several correlation relations have been found between the additional parameters, like the 
peak photon energy E Pt i, the variability V, the time lag ri ag , etc., and the GRBs radiated energy or luminosity. Assuming 
that GRBs emit radiation isotropically it is possible to relate L to the observed bolometric peak flux Ptoio and the bolometric 
fluence Sboio respectively by 

L = 4nd 2 L (z, cp)P bo i , (1) 
and 

E iso = 4TTd 2 L (z,c P )S boto (l + zy 1 , (2) 

where g(l is the luminosity distance, and cp denotes the set of cosmological parameters that specify the background cos- 
mological model. Equivalently one can use the total collimation corrected energy. It is clear from these relations that, in 
order to get the intrinsic luminosity, it is necessary to specify the fiducial cosmological model and its basic parameters. But 
we want to use the observed properties of GRBs to derive the cosmological parameters. Several procedures to overcome 
this complicated circular situation have been proposed (see for instance (Schaefer 2007; Basilakos & Pcrivolaropoulos 2008 
ICardone et al. 20081 IDem ianski et al. 2011)). In this paper we apply a Bayesian motivated technique, already implemented 
in ((pemianski ct a l. 2011[) to standardize 109 long GRBs with the well-known Amati relation, i.e. to find the E p ,j - _Bi so 
correlation parameters, where E P] i is the peak photon energy of the intrinsic spectrum and -Eiso is the isotropic equivalent 
radiated energy, in order to construct an estimated fiducial Hubble diagram that extends to redshifts z ~ 8, assuming that 
radiation propagates in a quintessential cosmological model. We begin our analysis using a minimally coupled self-interacting 
scalar field quintessense model, with an exponential potential. Parameters of this model are fixed by fitting appropriate es- 
timates to the set of type la supernovae data, the power spectrum of CMB temperature anisotropies and parameters of the 
observed large scale structure (see (|Dem ianski ct al. 2005)), which has been recently referred as fiducial model to standardize 
the GRBs and to construct the Hubble diagram using previous samples of GRBs compiled by Schaefer 2007, Dainotti et al. 
2008, Amati et al. 2008 - 2009, and Pemianski et al. 2011[) . Here we are extending this analysis by considering new updated 
dataset. Moreover we apply a local regression technique to estimate, in a model independent way, the distance modulus from 
the recently updated SNIa sample, referred to as Union (Amanullah ct al. 2010), containing 557 SNIa spanning the redshift 
range 0.015 < z < 1.55. The derived calibration parameters are used to construct a new calibrated GRBs Hubble diagram. 
We also compare the estimated and calibrated GRBs HDs. The scheme of the paper is as follows. In Section 2 we describe 
our statistical method to fit the E p ,i - Ei so correlation, to estimate the normalization and the slope of such a relation, and 
then construct the estimated and calibrated GRBs HD. In Section 3 we present an alternative procedure to calibrate the 
Amati relation in a cosmological-independent way, and we explore their implications for cosmography. Section 4 is devoted to 
discussion and conclusions. 



2 STANDARDIZING THE GRBS AND CONSTRUCTING THE HUBBLE DIAGRAM 

In this section we investigate the possibility of constructing the Hubble diagram from the E Pi i - Ei so correlation, here E Pi i is 
the peak photon energy of the intrinsic spectrum and Ej so the isotropic equivalent radiated energy. Ei so is defined by the Eq. 
([2]). This correlation was initially discovered in a small sample of BeppoSAX GRBs with known redshifts (Amati et al. 2002) 
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and confirmed afterwards by HETE-2 and SWIFT observations (|Lamb et al. 2005")) . fAmati 2006). Although it was the first 
correlation discovered for GRB observables it was never used for cosmology because of its significant "extrinsic" scatter. 
However, the recent increase in the efficiency of GRB discoveries combined with the fact that E Pi i - Ei so correlation needs 
only two parameters that are directly inferred from observations (this fact minimizes the effects of systematics and increases 
the number of GRBs that can be used by a factor ~ 3) makes this correlation an interesting tool for cosmology. Previous 
analyses of the E P} i - Ei ao plane of GRBs parameters showed that different classes of GRBs exhibit different behaviours, 
and while normal long GRBs and X-Ray Flashes (XRF, i.e. particularly soft bursts) follow the E p ^ - Ei ao correlation, short 
GRBs and the peculiar very near and sub-energetic GRBs do not (Amati et al. 2008). This fact may depend on the different 
emission mechanisms involved in different classes of GRBs and makes the E P: i - Ei BO relation a useful tool to distinguish 
between them (Antonclli 2009). The impact of selection and instrumental effects on the E P}1 - E- lso correlation of long GRBs 
was investigated since 2005, mainly based on the large sample of BATSE GRBs with unknown redshifts. Different authors 
came to different conclusions (see for instance (Ghirlan da at al. 20 05 )). In particular, (Ghirlanda at al. 2005) showed that 
BATSE events potentially follow the E P: i - Ei ao correlation and that the question to clarify is if, and how much, its measured 
dispersion is biased. There were also claims that a significant fraction of Swift GRBs is inconsistent with this correlation 
(Butler et al. 2007). However, when considering those Swift events with peak energy measured by broad-band instruments 
like, e.g., Konus-WIND or the Fermi/GUM or reported by the BAT team in their catalog (Sakamoto ct al. 2008) it is found 
that they are all consistent with the E Pt i - Ei so correlation as determined with previous /other instruments (1 Amati et al. 2 009). 
In addition, it turns out that the slope and normalization of the correlation based on the data sets provided by GRB detectors 
with different sensitivities and energy bands are very similar. These facts further support the reliability of the Amati correlation 
(Amati ct al. 2009). It is clear from Eq. © that, in order to get -Bi so , it is necessary to specify the fiducial cosmological model. 
In ([Dcmianski ct al. 201 lj , we fitted the Amati relation in a quintessence cosmological model, where the dark energy is 
described by the exponential potential of the scalar field discussed in (Dcmianski ct al. 2005). In our analysis we consider a 
sample of 109 long GRB/XRF, adding to the sample of 95 long GRB/XRF compiled in (Amati et al. 2008) and (Amati et al. 
2009) data of 14 unpublished GRBs, kindly provided by Amati in a private communication. Their redshift distribution covers 
a broad range of z, from 0.033 to 8.23, thus extending far beyond that of SNIa (z <~1.7), and including GRB 092304, the 
new high-z record holder of Gamma-ray bursts. 



2.1 Fitting the E P] i — Ei so relation and estimating its parameters 

In this section we present results of our analysis of the Amati correlation performed on a new updated dataset, assuming that 
the background cosmological model is one of the quintessence models that we have studied some time ago (Dcmianski c t al. 2005]) . 
and showed that it is consistent with the basic cosmological tests, and which has been recently referred as fiducial model to 
standardize the GRBs and to construct the Hubble diagram using previous samples of GRBs. First of all, we consider the 
E p i - Ei so relation in the form 



log ( ) = b + a log 

1 1 erg ' 



E p ,i(l + z) 



300 keV 



(3) 



where a and b are constants. In fitting this relation, we need to fit a data array {xi,yi} with uncertainties {<r x ,i, Vy,i}, to a 
straight line 

y — b + ax , (4) 

in order to determine the two fit parameters (a, b). Actually, the situation is not so simple since, both the (y, x) variables are 
affected by measurement uncertainties (a x ,a y ) which can not be neglected. Moreover, <J y /y ~ <j x /x so that it is impossible 
to choose as independent variable in the fit the one with the smallest relative error. Finally, the correlation we are fitting is 
not of theoretical nature, i.e, it is not (yet) derived from an underlying theoretical model determining the detailed features of 
the GRBs explosion and afterglow phenomenology. Indeed, we do expect a certain amount of intrinsic scatter, a int , around 
the best fit line that has to be taken into account and determined together with (a, 6) by the fitting procedure. Different 
statistical recipes are available to cope with these problems. As in (Dcmianski ct al. 2011), we apply a Bayesian motivated 
technique (D'Agostini 2005) maximizing the likelihood function C{a,b,cji„t) = exp [— L(a, b, <Ti nt )] with: 

L(a,b,a int ) = \ ^ In (of nf + a% i + a 2 ol i ) 
(yi - axi - b) 2 



(5) 



where the sum is over the ftf objects in the sample. Note that, actually, this maximization is performed in the two parameter 

d 

space (a, (Tint) since 6 may be estimated analytically by solving the equation —L(a, b, (Jint) = 0, as : 
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Figure 1. Marginalized likelihood functions: 
Ccr int , is obtained marginalizing over a. 



the likelihood function, C a is obtained marginalizing over <Xj„t; and the likelihood function, 
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Figure 2. Regions of 68%, 95% and 99% of confidence in the space of parameters a, Oi n t 



E; 



+ ait . + a <Ti 



+ a 2 ai 



(6) 



To quantitatively estimate the goodness of this fit we use the median and root mean square of the best fit residuals, defined 
as 8 — y bs — Vfit- To quantify the uncertainties of some fit parameter pi, we evaluate the marginalized likelihood d(pi) by 
integrating over the other parameter. The median value for the parameter pi is then found by solving : 



d(pi)dpi 



1 rPi 



Ci(pi)dpi 



The 68% (95%) confidence range 
Ci(pi)dpi = — — 



Pi,i,Pi,h) are then found by solving ( |D'Agostini 2 005) : 
Ci{pi)dpi , 



C t (pi)dpi 



Ci(p t )dpi 



(7) 



(8) 



(9) 



with e — 0.68 and e = 0.95 for the 68% and 95% confidence level. Just considering our correlation in Eq. ( [3| we find that 
the likelihood method gives a = 1.52, b = 52.67 and Oi nt = 0.41. In Fig. [2] we show the likelihood contours in the (a, (Tint) 
plane and in Fig. [3] we show the correlation between the observed log_E p ,i and derived logi?i so with our assumed background 
cosmological model. The solid line is the best fit obtained using the D'Agostini's method (D'Agostini 2005) and the dashed 
line is the best fit obtained by the weighted % 2 method. If one marginalize^ with respect to 6, then the likelihood values 
of a and &int are a — 1.52^q and (Tint — 0.45 i 0.05. The marginalized likelihood functions are shown in Fig. ^ The 
performed statistical analysis shows that the relation ([3]) has a statistical weight similar to the one exhibited by the other 
relations previously studied in (Dcmianski ct al . 2011[) . since both (5 me£ j and 5 rm s have almost the same values over the full 
set (of relations), where 8 — y i, s — Ufu- In our analysis we have assumed that the fit parameters do not change with the 
redshift, which indeed spans a quite large range (from z = 0.0331 up to z ~ 8). The limited number of GRBs prevents detailed 
exploration of the validity of this usually adopted working hypothesis, which we tested somewhat investigating if the residuals 



1 It is worth noting that in the marginalization procedure we have to take into account also the Eq. [6] 
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Figure 3. Best fit curves for the E Pi i — E- lao correlation relation superimposed on the data. The solid and dashed lines refer to the results 
obtained with the maximum likelihood and weighted x 2 estimator respectively. 
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Figure 4. Behaviour with the redshift (in a logarithmic scale, being u = log(l + z)) of the residuals, 5 = y b s ~ Vfit, f° r the E P}I 
Eiao correlation relation. We see that the figure does not exhibit any evidence of correlation. 



correlate with the resblft. We have not found any significant correlation, as shown in Fig. [4] Moreover we tested the fit of 
the E Pt i - Ei ao correlation relation with respect to the evolution with redshift, separating the GRB samples into four groups 
corresponding to the following redshift bins: z £ [0,1], z £ (1,2], z £ (2,3] and z > 3. We thus maximized the likelihood 
in each group of redshifts and determined the best fit parameters a, b at 65% confidence level, and the intrinsic dispersion 
Gint, as summarized in Table 1. It turns out that no statistical evidence of a dependence of the (a, 6, (Tint) parameters on the 
redshift exists. This is in agreement with what has recently been found by Ghirlanda et al. (2008) and Wang, Deng and Qiu 
(2008), or, as regards to other correlation relations, by Cardone et al. (2008) and Basilakos & Perivolaropoulos (2008), and 
also confirmed in our previous paper {Demianski et al. 2011[) . 



Table 1. Calibration parameters (a, b), intrinsic scatter Cj n t, median, 5 me£ j, and root mean square of the best fit residuals, <5 rms for the 
Amati correlation, evaluated in four ranges of redshift (z 6 (0,1], z £ (1, 2], z £ (2, 3] and z > 3, where however we have only 10 GRBs). 
Columns are as follows: 1. id of the reshift range; 2. maximum likelihood parameters; 3,4 68% confidence ranges for the parameters 
(o, (Ti n t)\ 5 ,6 median and root mean square of the residuals. 



Id 


(a, b,ai n t)ML 


(o_l CT , O+Ict) 


((fint) — Icti &int + ltT ) 




&rms 


z 6 [0, 1] 


(1.53, 52.6, 0.36) 


(1.28,1.69) 


(0.27,0.51) 


-0.07 


0.39 


z e (1,2] 


(1.38, 52.6, 0.55) 


(0.8, 1.9) 


(0.43,0.7) 


0.08 


0.58 


z e (2, 3] 


(1.38, 52.64, 0.56) 


(0.83,1.94) 


(0.43,0.75) 


0.08 


0.59 


z > 3 


(1.59, 52.9, 0.22) 


(1.08,2.01) 


(0.1 ,0.46) 


0.01 


0.29 
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Figure 5. The estimated quintessential Hubble diagram, with superimposed (solid line) the theoretical curve. 



2.2 Constructing the Hubble diagram 

Once the Amati correlation relation has been fitted, and its parameters have been estimated, we can now use them to construct 
the estimated fiducial GRBs Hubble diagram. Actually let us remind that the luminosity distance of a GRB with the redshift 
z may be computed as : 

"or (10) 

The uncertainty of dh{z) is then estimated through the propagation of the measurement errors on the involved quantities. In 
particular, remembering that our correlation relation can be written as a linear relation, as in Eq. (j4|, where y is the distance 
dependent quantity, while x is not, the error on the distance dependent quantity y is estimated as: 



trQogy) = \Ja 2 <r 2 (\ogx)+o'j nt , (11) 



and is then added in quadrature to the uncertainties of the other terms entering Eq. (|10|) to get the total uncertainty. The 
distance modulus fi(z) is easily obtained from its definition : 

fi(z) = 25 + 51ogd L (z), (12) 

with its uncertainty obtained again by error propagation. We finally estimate the distance modulus for each i-th GRB in 
our sample at redshift Zi, to build the Hubble diagram plotted in Fig. [5] In what follows we will refer to this data set as 
the fiducial GRBs Hubble diagram (hereafter, HD) since to compute the distances it relies on the calibration based on the 
fiducial quintessential model. We also investigated the impact of varying the parameters of our fiducial cosmological model, 
fitting the Amati correlation relation on a regular grid in the space of parameters of our quintessential model, as allowed by 
the data (Dcmianski ct al. 2005). For each point in the grid, we repeated all the steps described above to get the distance 
modulus to each GRB in the sample. We then collect these values and evaluate, for each GRB, the root mean square of the 
percentage deviation from the fiducial fi value. In this way we are performing a sort of average of the absolute percentage 
deviation, which allows us to quantify the order of magnitude of the studied effect. Actually it turns out that the distance 
modulus may be under or overestimated by a modest 0.3% with values never larger than 1%. 



2.2.1 Cosmological independent calibration: local regression on SNIa 

Although the above analysis has shown that the choice of the underlying cosmological model has only a modest impact on 
the final estimate of the distance modulus, we compared our estimated fiducial HD with a model independent calibrated 
HD, carried out using SNIa as distance indicators. We apply local regression to estimate the distance modulus fi(z) from 
the recently updated SNIa sample, the SCP Union2 compilation (Amanullah ct al. 2010), which is an update of the original 
Union compilation, now bringing together data for 719 SNe, drawn from 17 datasets. Of these, 557 SNe, spanning the range 
0.015 < z < 1.55, form the final sample considered in our analysis. To use this large dataset as input to the local regression 
estimate of fi(z) we have firstly to set a redshift Zi where /u(zj) has to be recovered and we order the SNIa dataset according to 
increasing value of \z — Zi\ and we select the first n = aAfsNia, where a is a user selected value and ATsnio, the total number 
of SNIa. Then we can fit a first order polynomial to the previously selected data, weighting each SNIa with the corresponding 
value of an appropriate weight function, like, for instance 

f (l-M 2 ) 2 M<i 

W{u) = I (13) 

[ o | u |>i, 
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Figure 6. Regions of 68%, 95% and 99% of confidence in the space of parameters a, <Ji„t, obtained by the local regression on SNela. 



Figure 7. Marginalized likelihood functions constructed by the local regression on SNela: the likelihood function, C a is obtained marginal- 
izing over ai n t; and the likelihood function, Co- int , is obtained marginalizing over a. 

and take the zeroth order term as the best estimate of fJ,(z). Here u — \z — Zi\/A and A is the maximum value of the \z — z,| 
over the subset chosen before. To estimate the error on we use the root mean square of the weighted residuals with 
respect to the best fit zeroth order term0. Having found an efficient way of estimating the distance modulus at redshift z in 
a model independent way, we can now fit the E P: i - Ei so correlation relation, using the local regression reconstructed fj,(z) in 
Eq. ©. We consider only GRBs with z < 1.55 in order to cover the same redshift range spanned by the SNIa data. For such 
subset of GRBs we apply the Bayesian fitting procedure described above to estimate the correlation parameters. We use the 
other GRBs to construct a new GRBs Hubble diagram that we call the calibrated GRBs HD. In Fig. [6] we show the likelihood 
contours in the (a,ai„t) plane, and the marginalized likelihood functions are shown in Fig. [7] 

Moreover it turns out that the estimated and calibrated HDs are fully statistically consistent, as shown in Fig. [8] 



3 AN ALTERNATIVE PROCEDURE TO CALIBRATE THE E PJ - Eiso RELATION IN A 
COSMOLOGICAL-INDEPENDENT WAY: IMPLICATIONS FOR COSMOGRAPHY 

In this section we propose a procedure to calibrate the E Pi i - Ei so relation using only the GRBs events with the redshift 
z < 1.55 without specifying a cosmological model. The luminosity distance estimations to GRBs are inferred from many 
known Type la supernovae, and based on an approximate formula for the luminosity distance which holds in any cosmological 
model, depending only on the shape of this function, more than on a power series expansion in the redshift parameter z 
(the coefficients of such an expansion being functions of the scale factor a(t) and its higher order derivatives), as in the 
cosmographic approach. Our starting point is the relation between the angular diameter distance Da and the luminosity 
distance Dl 

2 It is worth stressing that both the choice of the weight function and the order of the fitting polynomial are somewhat arbitrary. 
Similarly, the value of a to be used must not be too small in order to make up a statistically valuable sample, but also not too large 
to prevent the use of a low order polynomial. In our local regression routine we have performed an extensive set of simulations, ending 
up with a mock catalogue having the same redshift and error distribution of the actual SNIa one. This mock catalogue is used as input 
to the routine sketched above and finally the reconstructed fi(z) value for each point in the catalog are compared to the input one (see 
llCardone et al. 20081 for details). 
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Figure 8. Comparison of the distance modulus p.(z) for the calibrated and estimated GRBs Hubble diagram made up fitting the Amati 
correlation. The red solid line represent the graph of the function F(/isn) = PSN, and it turns out that it is also the fitting function if 
we do not include a constant term in the list of basis functions (nedeed to perform the fit procedure) . The dashed green line is the best 
fit function. 

D L = {l + zf Da, (14) 
where the angular diameter distance Da is a solution of the equation 
dz\ 2 d 2 D fd 2 z\dD 4tvG p 

with the following initial conditions: 
£>0)|*=o=0, 

(16) 

gLD(z) c 
dz Ho 

dx a 

In Eq. (1151) v is the affine parameter, T a g is the matter density tensor, k a = — — = — E, a is the vector field tangent to the light 

dv 

ray congruence, and E is the null surface along which the light rays propagate from the source. In the general form the Eq. (|15|l 
is very complicated. General properties of this equation have been extensively studied (see for instance, (Kantowskil998 
Kantowski, Kao & Thomas 2000; Kantowski & Thomas 2001; Dcmianski ct al.2003)). In most cases equation (|15[) does not 
have analytical solution, and from the mathematical point of view it can be reduced to a Fuchsian type with several regular 
singular points and a regular singular point at infinity. The solutions near each of these singular points, can be expanded 
in a series of hypergeometric functions. When we introduced the dimensionless angular diameter distance r = DHo/c we 
discovered (Dcmianski ct al.2003) that there is a simple function r(z), which quite accurately reproduces the exact numerical 
solutions of the equation fTS} for z up to very high values, it has the form 

r(z) = - (17) 



yfdi z 2 + (l + d 2 z + d i z 2 ) 2 



where di, d 2 and d$ are constants. Moreover the function (|17[) automatically satisfies the imposed initial conditions, so r(0) =0 
dv 

and = 1- This approximate expression immediately provides an empirical formula for the luminosity distance relation 

of the type la supernovae, through the Eq. {TJ. Fitting the relation (|17[1 to the SNIa Union dataset, we obtained the following 
best fit values for the fitting parameters di : 

di = -3.87 (18) 
d 2 = 1.51 (19) 
d 3 = 0.70. (20) 

In Fig.(]9j we show the approximate distance modulus with superimposed the SNIa data. The approximate function agrees 
with the real data within a relative error not larger than a few % in the SNIa redshift range, as shown in Fig. (llO|) . 

From the empirical approximate luminosity distance we can construct the empirical approximate distance modulus 
Ha PP rox{z), which we use to calibrate in a cosmological independent way, alternative to the one described above, the E p ^ 
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Figure 9. Behaviour of the approximate modulus of distance with the redshift, with superimposed the SNIa Union dataset. 
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Figure 10. Behaviour with the redshift of the relative residuals, 8 = — _J*t_^ Q f ^he gj. procedure which provides our empirical formula 
in Eq. dlTI . which agrees with the real data within a relative error of magnitude not more than few %. 



- Ei BO (Amati) correlation relation, following our standard procedure, and considering only the GRBs with z < 1.55 in order 
to cover the same redshift range spanned by the SNIa data. We use the other GRBs to construct a new GRBs Hubble diagram 
that we call the approximate calibrated GRBs HD. In Fig. [TT]we show the correlation between the observed log_E p ,j and de- 
rived log Eisc with our approximate luminosity distance. The solid line is the best fit obtained using the D'Agostini's method 
( |D'Ag ostini 2005]) and the dashed line is the best fit obtained by the weighted \ 2 method. In order to check the reliability 
of the approximate calibrated GRBs HD we compare it with the calibrated (SNIa) HD. It turns out that these HDs are fully 
statistically consistent, as shown in Fig. 1121 and 1131 and the resulting distances are strongly correlated with the Spearmans 
correlation p — 0.99. 
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Figure 11. Best fit curves for the E PiI — E- lso correlation relation, obtained using our empirical approximate luminosity distance formula, 
superimposed on the data. The solid and dashed lines refer to the results obtained with the Bayesian and Levenberg- Marquardt estimator 
respectively. 
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Figure 12. Comparison of the distance modulus fi(z) for the empirical approximate and calibrated (Snla) GRBs Hubble diagram. The 
datasct are fully statistically consistent are strongly correlated with the Spearmans correlation p = 0.99. The red solid line represent the 
graph of the function F(figjy) = p.SNi and it turns out that it is also the fitting function if we do not include a constant term in the list 
of basis functions (nedeed to perform the fit procedure). The dashed green line is the best fit function. 




3.1 Implications for high redshift cosmography 



Recently the cosmographic approach to cosmology gained increasing interest for catching as much information as possible 
directly from observations, retaining the minimal priors of isotropy and homogeneity and leaving aside any other assumptions. 
Actually, the only ingredient taken into account a priori in this approach is the FLRW line element obtained from kinematical 
requirements 



ds 



-c dt + a (t) 



dr 2 - 



kr 2 



+ r 2 dn 2 



(21) 



Using this metric, it is possible to express the luminosity distance ciz, as a power series in the redshift parameter z, the 
coefficients of the expansion being functions of the scale factor a(t) and its higher order derivatives. Such an expansions leads 
to a distance - redshift relation which only relies on the assumption of the Robertson - Walker metric thus being fully model 
independent since it does not depend on the particular form of the solution of cosmic evolution equations. To this aim, it is 
convenient to introduce the following parameters (Visscr 2004): 



H = if£. 

a dt 



(22) 
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* = (23) 

s = -f^ (25) 

(26) 

which are usually referred to as the Hubble, deceleration, jerk, and snap parameters, respectively. Their present day values 
(which we will denote with a subscript 0) may be used to characterize the evolutionary status of the Universe. For instance, 
go < denotes an accelerated expansion, while jo allows to discriminate among different accelerating models. It is worth noting 
that it is possible to infer implications for cosmography just using our empirical formula (|17[) of the luminosity distance in our 
analysis. Actually we first recast our approximate (Ll as a function of a new variable y = z/(l + z) (Vitagliano et al. 2010 
|Capo zzicllo a nd Izzo 2010[ ) in such a way that z £ (0, co) is mapped into y £ (0, 1), obtaining 

( \ 

(y) 



japprox 



C 

Ho~ 



(27) 



Expanding our approximate luminosity distance up to the fourth order in the {/-parameter, we get 
di ,c> . , . _ \ d 1 1 



d L (y) 



c 



y 



+ 4 - 4d 2 - da +6 +» 7;di(3d2 - 5) - d 2 + 5d 2 + 2d 2 (d 3 - 5) - 5(d 3 - 2) - (d 2 - 3)y + y } .(28) 



It is now possible to relate the di to the cosmographic parameters go, jo, so, by comparing the expansion in Eq. (|28[) with the 
standard expansion to the fourth order: 



d L (y) = — 

-no 



1 2 1 



3go 



jo] y + — [60 - 7jo- 



-10 - 32g + Wqojo + 6go + 21gg - 15gj] + s ] y 4 + 0(y 5 )} 
It turns out that 

1 + jo + go + 6 jo go - 2gp(l + 3go) + so 



dr = 
d 2 = 
d 3 = 



6(3 + go) 



3 + go 



14 + jo (5 - 4g ) + g (16 + 3(-l + g„)go) 



so 



(29) 

(30) 
(31) 
(32) 



12(3 + go) 

Inverting the systems of Eqs. (1301 1311 [32)) it is possible to recover the cosmographic parameters go, jo and so as functions of 
our fitted parameters di. Actually we get 

3 + 2d 2 (33) 

(34) 
(35) 



jo = 17 + 3di - 22d 2 + 6d 2 + 6d 3 

s = d 1 (51-24d 2 ) + 158d|-24di-8d 2 (35 + 9d 3 ) + 3(49 + 34d3). 



We actually get that 



These equations, together with the values of the fitting parameters in Eqs. f) 181 1191 120 p with the corresponding physically 
acceptable regions of confidenc^f] allow us to estimate the corresponding parameter confidence intervals for go, jo and so. 
/ go -1.11 0.27 \ 

jo —23.9 1.70 , the best fit being go = —0.45, jo = —12.3 and so = —99.3, which agree 
s -256.4 -34.7 J 

with the values found in literature (see for instance ( Vitagliano et al. 2010| |Capozzicllo and Izzo 2010)). In Fig. (|14[1 we 
plot our approximate cosmographic distance moduli together with the SNIa Union dataset. In addition, we investigate the 
possibility to use high redshift GRBs to determine parameters of our approximate cosmography. Therefore, in the following 



V 



3 Note that the use of the jerk parameter to discriminate between different models was also proposed in ( Salmi et al. 2003} in the context 
of the statefinder parametrization. 

4 It is worth noting that such physically acceptable regions are obtained imposing, on the confidence regions provided by the standard 
statistical fitting procedure, some priors, which have to guarantee that the approximate function in Eq. H17H preserves the special shape 
typical of the angular diameter distance. Instead, since the SNIa dataset used to fit the parameters di is limited in redshift by z < 1.39, 
the approximate function in Eq . I I 1 7B is not sampled at higher values of the redshift, and its behaviour could result misshaped. 
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we use the calibrated (with SNIa) GBRs HD, described above, to fit the di parameters and then to derive the cosmographic 

/ di 0.950028 11.2294 

parameters go, jo, and so- We obtain the following (2a) parameter confidence intervals 



and 



d 2 -2.76299 1.92514 
y d 3 0.140305 0.836675 

it turns out that the corresponding parameter confidence intervals for go, Jo and so (through the Eqs. 1301 1311 I32[) are: 
go -1-09 0.10 \ 

jo —0.3 2.71 , the best fit being go = —0.51, jo = 0.69 and so = —0.48, in agreement with other results 
s -1.23716 8.34 J 

in literature (see for instance ( |Vitagliano et al. 2010| IGao et al. 2010)1 ). In Fig. (|15[> we plot our approximate cosmographic 
distance modulus together with the SNIa Union dataset and the calibrated GBRs HD. The reliability of the reconstruction 
is measured by the relative residuals e = — Uhl^ shown in Fig. [TB] We finally perform our approximate cosmographic 
analysis, considering a whole dataset containing both the SNIa Union dataset and the calibrated GBRs HD, which we call the 
cosmographic dataset. 




Figure 16. Residuals (relative) e : 



for our approximate cosmography: the empty black boxes correspond to the Snela Union 



dataset and the filled red diamond to the calibrated GBRs HD. 
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Figure 17. Distance modulus for the best-fit values of our cosmography performed with the cosmographic dataset. The filled black 
circles correspond to the SNIa Union dataset and the empty red triangles to the calibrated GBRs HD. 



50 - 
45 - 

a. 

40 : 
35 : , 




0.0 0.2 0.4 0.6 0.8 



y 

Figure 18. Distance modulus for the best-fit values of the standard cosmography performed with the cosmographic dataset. The filled 
black circles correspond to the SNIa Union dataset and the empty red triangles to the calibrated GBRs HD. 



We obtain the following parameter confidence ( at 2a) intervals for the d; and the cosmographic parameters respectively: 
I d 1 -1.35 2.26 \ 

d 2 0.71 1.30 , 
v d 3 0.39 0.59 J 

I q -1.58 -0.39 \ 

jo 2.67 8.87 . The best fit being qo = —0.9, jo = 5.25 and so = 27.25, which again agree with the values 
v s -6.88 46.9376 J 

found in literature. In Fig. [17] we show the approximate cosmographic distance modulus together with the cosmographic 
dataset. Moreover we observe that since the cosmographic dataset includes also the GBRs HD, which spans a quite large 
range of redshift up to z ~ 8, the approximate function in Eq.{T7]) is sampled also at higher values of the redshift, and its 
behaviour is not misshaped, even without any prior on the confidence regions, as nedeed above. In order to further check 
the reliability of our procedure, we compare the results summarized above with that provided by the standard cosmography. 
We apply this procedure to the cosmographic dataset only, and obtain the following parameter confidence (2a) intervals: 
( q -1.09 -0.08 \ 

jo —0.30 2.71 , which fully agree with the previous results. Also from the point of view of the analysis of residuals 
v s -1.24 8.34 J 

it turns out that our approximate cosmography is statistically fully consistent with the standard cosmography: actually 
we obtain a fully compatible values for the rootmean square and the correlations in both cases. In Fig. [18] we show the 
standard cosmographic distance modulus together with the cosmographic dataset. In short, let us note that we implemented 
a cosmographic analysis, using both the supernovae and the Gamma Ray Bursts data, which allow us to obtain constraints 
on the parameters of cosmography (starting from our approximated luminosity distance), which we also compare with the 
results of a standard approach. Therefore our approach is very different from the one used in (Capozzicllo and Izzo 2010), 
where they obtain a cosmographic luminosity distance in the y-redshift, which is used to calibrate the E Pi i - Ei so relation 
using a weighted x 2 estimator, but without making up the GRBs Hubble diagram. Moreover in that case the constraints on 
the cosmographic parameters are obtained from the SNIa data only. 
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4 DISCUSSION AND CONCLUSIONS 

Recently several interesting correlations among the Gamma Ray Burst (GRB) observables have been identified. Proper 
evaluation and calibration of these correlations are needed to use the GRBs as standard candles to constrain the expansion 
history of the universe up to redshifts of z ~ 8. Here we used the GRB data set recently compiled by Amati et al. (2008 
- 2009) to investigate, in a quintessential cosmological scenario, the E Pi i - E{ ao correlation relation. Marginalizing over the 
normalization, b, our Bayesian analysis provides the following parameter confidence (at 3ct) intervals for the the best-fit power- 



coefficient turns out to be b = 52.67. We used the fitted parameters to construct an estimated fiducial Hubble diagram 
that extends up to redshifts z ~ 8. In the first part of our analysis we have assumed that the fit parameters do not change 
with the redshift. The limited number of GRBs prevents detailed exploration of the validity of this usually adopted working 
hypothesis, which we tested somewhat investigating if the residuals correlate with the redshift. We have not found any 
significant correlation. Moreover we tested the fit of the E p ,i - E{ ao correlation parameters with respect to the evolution with 
redshift, binning the GRBs into four groups with redshift from low to high, each group containing a reasonable number of 
GRBs. We thus maximized the likelihood in each group of redshifts and determined the best fit calibration parameters. Our 
analysis indicates that the presently available GRBs datasets do not show statistically unambiguous evolutionary effect with 
the cosmological redshift, but this should be tested further with larger GRBs samples. We also investigated the impact of 
varying the parameters of our fiducial cosmological model, fitting the Amati correlation relation on a regular grid in the space 
of parameters of our quintessential model, and evaluating, for each GRB, the root mean square of the percentage deviation 
from the fiducial fj, value. In this way we performed a sort of average of the absolute percentage deviation which provides 
the variation of the fitting parameters a and b. It turns out that the distance modulus may be under or overestimated by a 
modest 0.3% with values never larger than 1%, so that the underlying cosmological model (in terms of varying the values 
of its characteristic parameters) has only a modest impact on the final estimate of the distance modulus. However, this fact 
does not imply that the cosmological constraints that can be obtained from GRBs data with this method are also marginal: 
already in IjDem ianski et al. 201 Q we tested the stability of values of the fitted correlation parameters by considering an 
ad hoc definition of the luminosity distance that gives much larger distances to objects at z > 2 than either of the fiducial 
model. It turned out that such artificial crazy luminosity distance is changing the values of the correlation parameters, but 
the difference is not dramatic at the la confidence level. However the situation changes when we consider, in the space of 
parameters, the regions of la, 2 a and 3<r of confidence for our crazy model. It turns out that with respect to the same 
regions constructed for the fiducial model they overlap only at 1 a, but differ consistently at higher levels of confidence. As 
a consequence, when we made up the GRBs Hubble diagram the number of GRBs deviating from the fiducial \x more than 
2a increases from N = 14, in the case of the quintessential cosmological model, to N = 51 in the case of the crazy model 0. 
The calibration of the Amati correlation, as well as the other known correlations, in several cosmological scenarios is therefore 
already needed to fully use the GRBs as a cosmological probe. Actually, quite recently, in (Dia ferio et al. 2011[) the GRBs are 
used to test the ACDM vs. conformal gravity; and in a forthcoming paper we are standardizing the GRBs with the Amati 
relation to test cosmological models based on extended theories of gravity. Moreover we apply a local regression technique to 
estimate, in a model independent way, the distance modulus from the recently updated SNIa sample. The derived calibration 
parameters are used to construct an updated GRBs Hubble diagram, which we call the calibrated GRBs HD. We also compare 
the fiducial and calibrated GRBs HDs, which turned out to be fully statistically consistent, thus indicating that they are 
not affected by any systematic bias induced by the different calibration procedures. This means that the high redshift GRBs 
can be used to test different models of dark energy settling the circularity problem. Finally we propose a new approach to 
calibrate the GRB relations in a cosmologically independent way, by using an approximate luminosity distance relation, which 
holds in any cosmological model. We use this calibration of the Amati relation to construct an empirical approximate HD, 
which is fully consistent with the calibrated GRBs HD, probing in such a way the reliability of our approach, which could 
provide a robust procedure to calibrate in a cosmological independent way different GRBs correlation relations, specially if 
the available dataset are poorly populated in the SNIa range of redshift, so that, in such a case, the SNIa fitting procedure 
fails. We finally investigated the implications of such an approach for the high redshift cosmography. Actually, starting from 
the estimation of the constant di , di and dg present in our approximate luminosity distance relation, we constructed the map 
which connects our dts to the parameters describing the kinematical state of the universe qo, jerk jo, and snap so. This map 
is a core for our approximate cosmography , which we actually applied to a whole dataset containing both the SNIa Union 
dataset and the calibrated GBRs HD, which we call the cosmographic dataset. Using the cosmographic dataset, we show that 
the deceleration parameter qo up to the 3a confidence level is definitively negative. The constraints for the jerk and snap 




1.263 1.77 
0.36 0.5, 




5 In (Dcmianski ct al. 2011 ) the results concerning the crazy model are inferred from fitting five correlation relations different from the 
Amati relation. However we tested that for the Amati relation also we obtain the same kind of behaviour. 
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parameters are, instead, less strong. Finally it turns out that our results are fully consistent with the results obtained using 
the standard approach, showing that our new approximate cosmography allows robust results in a wide range of redshifts. 
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